Experimental study of the jamming transition at zero temperature 
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We experimentally investigate jamming in a quasi- two-dimensional granular system of automati- 
cally swelling particles and show that a maximum in the height of the first peak of the pair correlation 
function is a structural signature of the jamming transition at zero temperature. The same signature 
is also found in the second peak of the pair correlation function, but not in the third peak, reflecting 
the underlying singularity of jamming transition. We also study the development of clusters in this 
system. A static length scale extracted from the cluster structure reaches the size of the system 
when the system approaches the jamming point. Finally, we show that in a highly inhomogeneous 
system, friction causes the system to jam in series of steps. In this case, jamming may be obtained 
through successive buckling of force chains. 

PACS numbers: 81.05.Rm, 61.43.Fs, 83.80.Fg 
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I. INTRODUCTION 

Glasses return to the liquid state upon heating - they 
become soft and can flow. Sand flowing through a pipe 
or out of an orifice can easily jam and become rigid. In 
fact, large classes of materials, ranging from polymer 
melts to foams, from glasses to dense colloidal suspen- 
sions and granular matter, show a similar transition be- 
tween a flowing liquid-like state and a non-equilibrium, 
disordered soHd state. How to understand the nature 
of this transition is one of the central questions in dif- 
ferent fields of materials science [H-lsl. The recently 
proposed jamming phase diagram provides an approach 
to unify this effort [gI, It^. In such a phase diagram 
(Fig. [Tj), athermal systems such as granular media sit in 
the inverse-density (1/0) /shear- stress (S) plane, and ther- 
mal systems, such as glass- forming liquids, sit in the 
inverse-density (1/0) /temperature (T) plane. When de- 
creasing the temperature, decreasing the shear stress or 
increasing the density, a system goes from an unjammed 
phase into a jammed state. At the jamming transition, 
a material becomes rigid and loses its ability to explore 
the entire phase space efficiently; it therefore falls out of 
equilibrium. 

Due to its clean and clear signature, jamming along the 
axis in the jamming phase diagram, at zero temper- 
ature and zero shear stress, has so far attracted the most 
attention in simulations and theoretical studies 
At the zero-temperature jamming point for frictionless 
spheres (called point J) one observes the onset of rigid- 
ity along with other interesting phenomena such as the 
divergence of the pair correlation function and the ap- 
pearance of anomalous soft modes (7l-[l3[ flBj . It has been 
hypothesized that the properties of point J will influence 
the jamming transition also nearby, similar to critical 
points in second order phase transitions. Indeed, it has 
been found recently that the signature of the jamming 
transition at zero temperature shows its vestige at finite 
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FIG. 1: (Golor online) Jamming phase diagram 0,01 • The 
volume bounded by the black lines near the origin is the 
jammed phase. When a material crosses the black line from 
the outside, it goes through the jamming transition. The jam- 
ming point along the 1/0 axis, point J, is marked by a red 
dot. 



temperature [isl, |l9| . Although there has been intensive 
theoretical and simulation work, only a few experiments 
have been conducted to study the nature of the zero- 
temperature jamming transition (2qI-[22|. Especially, the 
structural signature of the jamming transition has not 
been directly addressed in experiments by far. Corwin et 
al. probed the structural signature of the jamming tran- 
sition indirectly at zero temperature along the S axis of 
the jamming phase diagram [23[. Can one directly ob- 
serve the structural signature of a zero-temperature jam- 
ming transition along the 1/0 axis experimentally? How 
is the picture of jamming modified in a real system with 
frictional interactions? This paper investigates these es- 
sential questions. 

Here we study the zero-temperature jamming transi- 
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tion experimentally in a quasi-two-dimensional granular 
system of macroscopic particles. By continuously and 
uniformly increasing the packing fraction, the system 
is forced to go through the zero-temperature jamming 
point. In this process, we find a structural signature, 
which is shown as a maximum of the height of the first 
peaks of pair correlation function. A similar signature 
was previously seen in a colloidal sample at non-zero ef- 
fective temperature isl, [l9|. Here, we observe a maxi- 
mum in the second peak as well and, by measuring the 
pressure on the boundary of system, we show that these 
features coincide with the onset of rigidity. We also find 
that en route to the jammed phase, our athermal system 
always self-organizes itself into a structure consisting of 
particle clusters. A static length scale can be extracted 
from this structure; it varies between a few particles to 
the size of the system when the jamming point is ap- 
proached. Finally, we show that friction can result in 
multiple jamming points in the presence of highly het- 
erogeneous particle arrangements. A system with friction 
shows a clear historical dependence. Frictional effects are 
shown to be reduced or eliminated when small amplitude 
vibrations are introduced. 



II. EXPERIMENTAL METHOD 

The granular material used consists of tapioca pearls, 
which are spherical particles made of starch (Fig. in- 
set). The average diameter of dry pearls is 3.3mm and 
the cumulative distribution of the sizes of dry pearls is 
shown in Fig. [2^. The polydispersity of the pearls is 8%. 
Such a spread in size is essential to avoid crystallization in 
a two-dimensional (2D) system. One important property 
of tapioca pearls is that, when submerged in water, they 
uniformly expand in size (Fig.[2]3). The final diameter of 
a tapioca pearl can be 1.7 times larger than its original 
value. The swelling process is very slow. It takes 24 hours 
for a pearl to reach its final, fully swelled state. Thus, 
the system is quasi-static. During swelling, particles keep 
their approximately spherical shape. The contact inter- 
action between fully swelled particles is purely repulsive 
and of Hertzian type (Fig. [2t). The strain versus stress 
curve on an individual particle is measured with an In- 
stron System (Model 5869): a particle is put between two 
horizontal metal plates and the normal force on the top 
plate is measured while the gap between the two plates 
decreases. For particles which are not fully swelled, the 
inter-particle contact is still purely repulsive but may de- 
viate from Hertzian form. Some particles may also show 
slightly plastic deformation. 

At the beginning of an experimental run, we put tapi- 
oca pearls randomly into a square cell submerged in wa- 
ter. The schematics of the setup can be seen in Fig. [H 
Two different sizes of cells are used. The side length of 
the larger cell is L = 54.6 cm, which can hold over 15,000 
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FIG. 2: Properties of tapioca pearls, (a) Cumulative distri- 
bution of diameter of dry tapioca pearls. The vertical dashed 
line marks the average diameter of pearls, (b) Diameter of 
a swelling pearl under water relative to its original size as a 
function of time. Inset: optical image of dry tapioca pearls. 
The white scale bar is 1 cm. (c) Compression force on a fully- 
swelled tapioca pearl as a function of its strain. The strain is 
defined as Ax/d, where Ax is the deformation of the particle 
under compression and = 4.7 mm is the original size of the 
swelled particle at zero compression. The slope of the dashed 
line is 3/2. 
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FIG. 3: (Color online) Schematics of the experimental setup. 
Shown here is the cross section of the cell. The blue (gray) 
area indicates water. 

particles with initial packing fractions (/)initiai ^ 0.63 
(among which about 10,000 particles are studied in the 
central area to avoid the boundary effect); the side length 
of the smaller cell, which can contain roughly 2,000 par- 
ticles with about 1,000 particles in the central area, is 
L = 16.5 cm. In the small cell, we installed a force sen- 
sor (Futek load cell. Model LSB200) along one of its sides 
so that the force or pressure along that boundary can be 
measured. For both the large and small cells, the gap 
between the top and bottom plates is kept by a 5.2 mm 
spacer enclosing the entire boundary. To allow water to 
flow in and out of the cell, a few thin washers (0.254 mm 
in thickness) are put on top of the spacer. The total gap 
thickness (spacer + washer = 5.454 mm) prevents fully 
swelled particles from buckling out of the plane signifl- 
cantly to form two layers. The cell can also be coupled 
to a mechanical shaker at the bottom, which can vibrate 
the entire cell vertically. However, unless explicitly indi- 
cated, the experiments described here are done without 
shaking. Vibration is only introduced at the end of Sec. 
V for controlling the friction between particles. An im- 
age of the particles is recorded by a camera mounted 
above the cell. For a system to reach its fully jammed 
stationary state, a typical experiment takes between 17 
and 24 hours. An image was taken every 20 s so that 
3000 to 4000 images were recorded for each experiment. 
From these images, the center and trajectory of individ- 
ual particles were extracted. 

We track the center of particles based on an algorithm 
developed for colloidal systems [24] . The images obtained 
are flrst processed with a band-pass fllter to ehminate any 
global intensity variation and the pixel-size noise. This is 
a standard procedure before tracking the center of parti- 
cles [24]. The resulting images are then convoluted with 
a mask consisting of a white disk with radius a little 
smaller than the size of particles. This eliminates any 
intensity variation on the top surface of individual parti- 



cles. The local maxima are then located. The center of 
particles can be found more precisely by calculating the 
centroid of a blob around each local maximum . The 
center of particles found in this way is more accurate at 
early times or low packing fractions. When the system 
is deep inside jammed phase, the interface between par- 
ticles has much lower contrast. Therefore, the error of 
particle tracking becomes larger. When (j) > 0.90, about 
1 — 2% of particles are missed by the algorithm. 

There are several advantages of this system. First, dif- 
ferent from other granular systems for studying the jam- 
ming transition, where the packing fraction is changed by 
either changing the number of particles or by changing 
the volume of system from the boundary [20*, 25-^7], here 
we can continuously and uniformly increase the packing 
fraction across the entire system. Second, the system is 
quasi-static due to the extremely slow swelling of the par- 
ticles. This allows the static structure of the pack to be 
easily investigated. Third, after swelling tapioca pearls 
are much softer than other commonly used granular ma- 
terials such as glass beads. By assuming the Poisson's ra- 
tio of particles around 1/3, we can estimate the Young's 
modulus of swelled particles from the force-strain curve 
shown in Fig. [2t. The Young's modulus of swelled tapi- 
oca pearls is 0.060 ± 0.005 GPa, which is three orders of 
magnitude smaller than that of glass beads. Thus the 
system can reach far inside the jammed phase. This is 
essential to directly see any small structural signatures 
of the jamming transition. With hard granular materi- 
als, structural signatures can only be probed indirectly 
by measuring the contact force distribution between par- 
ticles [235. Also, since the entire system is under water, 
the friction between particles and the bottom plate of 
the cell is reduced due to lubrication. It is interesting to 
note that our system is a 2D version of the old experi- 
ment done by Stephen Hales in 1727 [28, 29]. In order 
to flnd out how many contacting neighbors a spherical 
particle has in a dense pack. Hales put peas into a flxed 
volume container full of water and counted how many 
dimples each pea had after they had swelled. 



III. STRUCTURAL SIGNATURE OF JAMMING 
TRANSITION AT T = 

We initially prepared the 2D samples in a dense, but 
unjammed phase (Fig. |4^). As the particles become 
larger, the packing fraction, 0, of the pack increases uni- 
formly across the entire system (Fig. |4)3 andllJ:). At a 
certain moment the system crosses the jamming point 
and goes into the jammed phase (Fig. |4]i). The jamming 
transition appears to be continuous. The question posed 
here is whether one can identify the jamming point by 
looking merely at the structure of the pack. 
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FIG. 4: Jamming transition at zero temperature. The di- 
ameter of the particles, and therefore the packing fraction, 
increases with time, (a) (j) = 0.62 at t = h, (b) = 0.76 at 
t = 1.10 h, (c) (j) = 0.84 at t = 2.81 h, and (d) (p = 0.92 at 
t = 16.66 h. The white scale bar is 2 cm. 



A. Experimental results 

To study the structure of the pack, we measure its 
pair correlation function, g{r) [29|. As shown in Fig. [5l 
at each ^, g{r) has an oscillating shape characteristic of 
any disordered medium. As (f) increases from the un- 
jammed phase, the height of the first peak of g{r), ^i, in- 
creases first. However, when (j) is above (j)c = 0.84 ±0.02, 
gi begins to decrease (Fig. [6j). Here, we measured the 
packing fractions from the two-dimensional projection of 
images. Hence, the average size of particles at (j)c can be 
estimated as d = <io(0c/0initiai)"'^^^, where do = 3.3mm 
is the average initial size of particles. We identify (j)c by 
fitting ^i(0) near its peak with a peak function. Thanks 
to symmetric shape of ^i(^) near a Gaussian func- 
tion provides a good fitting. Presumably any other sim- 
ilar peak functions would lead to the same value of ^c- 
The non-monotonic trend in gi (0) indicates a structural 
signature. Is this signature a signature of jamming tran- 
sition? In other words, is 0c the jamming point? 

To answer that, we checked two well-defined criteria 
for jamming. First, the jamming point is supposed to 
mark the onset of rigidity in a disordered system [7]: a 
system composed of particles with finite range purely re- 
pulsive interactions at T = begins to build up pressure 
on its boundary at this point. As shown in Fig. [7^, the 
force measured at the boundary of the cell is zero ini- 
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FIG. 5: (Color online) Pair correlation, g{r), at different pack- 
ing fractions. The distance r is given in the unit of the average 
dry particle diameter, do. The height of the flrst peak of g{r), 
gi, is indicated. 
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FIG. 6: Structural signature of the jamming transition. The 
height of the first peak of the pair correlation, ^i, is plotted 
as a function of 0. The position of the maximum is indicated 

as (be- 



tially when the system is unjammed, and begins to de- 
viate from zero at a ^ consistent with the peak of g\. 
To illustrate the detail of the onset of the jamming, we 
also plot the force in the logarithmic scale (Fig. [Jt)). Be- 
low the jamming point, the force fluctuates around the 
noise level (below 100/iN) of the instrument. As ^ in- 
creases further, the slope of the curve changes sharply 
at ^ = , which is indicated by the left red arrow in 
Fig. [7)3. Above 0f , the pressure on the boundary in- 
creases signiflcantly. Note that 0f is a little ahead of 
but very close to ^c- The phenomenon is robust for ex- 
periments with uniform initial packing fractions. We flnd 
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that (l^c/^'c — 1-009 ± 0.006, which clearly suggests that, 
0c 5 and therefore the structural signature we found in 
is directly related to the jamming transition. We 
suggest that the small difference between (j)^ and (j)c is 
due to friction in the system. Particles compressed onto 
the force sensor can be held in a force balance by friction 
with other particles and with the bottom of the cell, and 
therefore are not jammed globally with all the particles 
in the system, which always results in a smaller cf)^ than 
0c- However, the static friction in the aqueous system 
near the isostatic point of the jamming transition is too 
small to sustain much stress from the swelling of parti- 
cles. The force balance maintained by friction will break 
down quickly. Therefore, 0c — 0f 1- 

Another supporting evidence is from the motion of par- 
ticles. Even though the system is athermal, a particle can 
still be displaced when it touches other particles during 
the swelling process (Fig. [8^ andlSb). We shall discuss 
the displacement of particles in more detail in Sec. IV. 
Here, it is sufficient to know that the motion of particles 
stops at the jamming point, which is self-evident as the 
kinematic criterion of the jamming transition. Ideally, 
in a homogeneous system, swelling particles will touch 
the boundaries of the cell in different directions simulta- 
neously, which results in an extremely sharp drop of the 
average velocity of particles at = 0c- However, in an in- 
homogeneous frictional system, some particles may reach 
a boundary of the cell faster and stop the motion ffist, 
while particles in other parts of the system still move. 
Therefore, there exists a finite interval A0 for diminish- 
ing of particle motion. For an experiment with a uni- 
form initial condition, the interval A0 is small. To show 
the average behavior of particles' motion quantitatively, 
we measured the mean square displacement of particles, 
(1)^(0)) = ^iLi experiments. Here, 

Di{t{(j))) = fi{t -\- At) — fi{t) is the displacement of the 
particle i at time t within a small time interval At (or 
equivalently at packing fraction within an interval A0 
since the monotonic dependence of t(0)) and the summa- 
tion is over all the particles in the studied area. fi{t) is 
the center of the particle z at We fixed At to be a small 
constant compared with the time scale of the swelling of 
particles, so Di(t)/ At is approximately the instantaneous 
velocity of particle z at t. Hence, the mean squared veloc- 
ity should show qualitatively the same behavior as that 
of {D'^). Note that the mean square displacement defined 
here is different from the more common concept used in 
studying the diffusion of Brownian particles. The result 
of (1)^(0)) is shown in Fig. [Hb. As one can see from 
the plot, (1)^(0)) decreases to zero quickly in an inter- 
val between 0^ and 0r as indicated in the figure. For 
all experiments with uniform initial packing fraction, we 
found that A0/0c = (0r - 00/^c = 0.03 ± 0.01. More 
importantly, 0c is always located between 0^ and 0^ (i.e. 
<j)i < (j)c < (t^r)- This conffims the argument that the 
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FIG. 7: (Color online) Comparison of the structural signature 
5'i(0) with the force measured along the boundary F(0) (a) in 
a linear- linear plot and (b) in a logarithmic-linear plot. The 
black squares show 5'i(0) on the left and the red line shows 
F(0) on the right. The vertical dashed line marks 0c and the 
horizontal dashed line indicates zero force. The onset of force 
0f is indicated in (b). 

structural signature in ^i(0) is due to the jamming tran- 
sition. 

Both the mechanical and kinematic measurements 
show that 0c indeed corresponds to the jamming point of 
frictional systems. Therefore, the peak in ^i(0) manifests 
the structural signature of the zero-temperature jamming 
transition. But how can we understand this structural 
signature? What is the physical or geometric origin of 
it? We shall discuss it next. 

B. Discussion 

It is well known that when a molecular liquid goes 
through glass transition upon decreasing T in the jam- 
ming phase diagram, no structural signature can be ob- 
served [2I, [32I; gi monotonically increases as the sys- 
tem passes through the glass transition temperature Tg. 
Then why does gi show a peak when the system goes 
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FIG. 8: (Color online) Comparison of the structural signa- 
ture with the mean square displacement {D'^{(j))). (a) 
and (b) are sketches to show particle motion. If two particles 
touch (solid circle), the centers of particles will be displaced 
after they swelled up (dashed circle), (b) Particles sitting at 
the edge of a cluster will be displaced most. The displace- 
ment is proportional to the size of cluster, (c) Comparison 
between gi{(j)) and (D'^ {(/))). The displacement of particles D 
is measured in a time interval At = 800 s. The black squares 
show gi{(l)) on the left and the red circles show (D^^cj))) on 
the right. The interval of decreasing {D'^ ((/))) is indicated by 
(f)i on the left and by (f)r on the right. The vertical dashed line 
indicates 0c. Due to random noise in the system, (D^((/))) 
does not approach zero at large 0. The source of the noise is 
discussed in the main text of Sec. IV. 



through jamming transition along the l/cj) axis? This 
was discussed in Ref. ^^8], where a vestige of the zero- 
temperature signature was observed in a finite temper- 
ature colloidal system. When (j) approaches (j)c from 
the unjammed side, particles are pushed closer to each 
other. Therefore, ^i, which indicates the probability that 
nearest-neighbors of a particle are located at the same 
distance, increases as the total number of nearby neigh- 
bors increases. However, above the jamming transition, 
particles begin to overlap and deform. Since the degree 
of deformation depends on the local environment of a 



particle, the distribution of distances between two parti- 
cles in contact becomes broader. Meanwhile, the number 
of nearest-neighbor particles (the coordination number) 
does not increase appreciably, i.e., the area under the 
first peak of g{r) remains roughly constant. As a result, 
^1, the height of the first peak in ^(r), decreases. In 
simulations with monodisperse frictionless particles, at 
the jamming point, all particles are precisely one particle 
diameter away from their nearest neighbors. Hence, gi 
diverges at this point [9]. In our system, the polydisper- 
sity of the particles reduces gi. Even without friction, at 
the jamming point, the distribution of distance between 
two particles in contact still has a finite width reflecting 
the size distribution. 

Friction may complicate the situation further. Both 
the coordination number and the deformation of parti- 
cles under compression are profoundly changed in the 
presence of friction as indicated in previous simulation 
and theoretical works 
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Hence, it is not 
straightforward to extend the results of ideal frictionless 
system on which most simulation and theoretical stud- 
ies focus [H, [13, Esl, [13 to experiments with real 
frictional granular matter. Experimentally, Majmudar et 
al found that the increasing of the coordination number, 
Z — Zc^ and the pressure, P, of the system as a function 
of ^ — agrees with the mean-field theory for friction- 
less particles [2o|. However, the two experiments on the 
sound propagation near the surface of loosely compacted 
granular packs show that the ratio of the shear modulus 
to the bulk modulus, G/B^ stays constant rather than di- 
minishing as the pressure of the system approaches zero 
[21I . [22! . which contradicts to the result of frictionless 
particles [3, [l3- Until now, no direct measurement has 
been conducted on the structural signature of the jam- 
ming transition in a real granular system. Here, we show 
that structural signature predicted with frictionless par- 
ticles [9] persists in the system of real granular matter 
with frictional contact, although the signature is modi- 
fied significantly. 

Finally, it is also interesting to look for this structural 
signature in the other peaks of the pair correlation func- 
tion. The system with the larger cell contains enough 
particles to show the first four peaks of g{r) clearly 
(Fig. [5]). As seen in Fig.[9K, the height of the second peak 
of ^(r), ^2, as a function of (j) also shows a maximum at 
the same value of 0c , at which gi{(j)) shows a maximum. 
However, the amplitude of this maximum is smaller than 
that of gi (Fig. [Hh). One may think that the height of 
the nth peak of the pair correlation function, ^f^, would 
also show a peak at 0c but with decreasing amplitude as 
n increases. However, when we measured the height of 
the third peak, ^3(0), no peak is found: instead gs in- 
creases rapidly at small and slows down or plateaus in 
some cases at large (Fig. [9)3). The structural signature 
of jamming transition apparently manifests itself in the 
first and the second peaks of the pair correlation function 
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but not in the peaks at larger separation. As a compar- 
ison, in the simulations with frictionless monodisperse 
particles, the second peak of the pair correlation func- 
tion splits into two sub-peaks [9|, |29j both of which have 
a divergent slope ^9|]. Furthermore, no singular behavior 
of gs is found in simulations So far, unlike the case 
of gi , no clear geometric picture exists for why g2 should 
diverge at the jamming point [9]. It is hypothesized that 
the underlying mechanism for the singularity of simula- 
tions may be related to the experimental finding. 



IV. FORMATION OF STATIC CLUSTER 
STRUCTURE 



In this section, we investigate the displacement field of 
the particles. We shall show that, for any typical pack, 
the system eventually organizes itself into clusters. One 
can extract a static length scale from the sizes of the 
clusters. This length grows dramatically from the size of 
a few particles to the size of the entire system when the 
system approaches the jamming point. 



A. Displacement of particles 

The centers of two contact particles separate due to the 
enlargement of their radiuses as illustrated in Fig. [8^. 
Thus, by tracking centers of particles, one can observe 
clear motion in this athermal system. 

Now let us first have a look at the average behavior of 
particles at different stages through the jamming tran- 
sition. As shown in Fig. [8h, the mean square displace- 
ment of particles, (1)^(0)), shows a non-monotonic be- 
havior. Since only a few pairs of particles have contacts 
initially, {D'^{(j))) is small at the beginning (Fig. [St). As 
the size of particles increases, more particles form con- 
tacts and are displaced; {D'^{(p)) increases correspond- 
ingly. When the system approaches the jamming point, 
almost all the particles join into a contact network so 
that {D'^{(j))) reaches a maximum. However, shortly af- 
ter that, the moving particles begin to touch the bound- 
ary and {D'^ {(/))) quickly drops toward zero as the entire 
system jams. The position of the peak {D'^{(f))) is always 
before the peak of ^i(^). Deep inside the jammed regime, 
there are rare buckling events: particles suddenly change 
their relative positions on a much shorter time scale than 
that for particle swelling. These events are discrete and 
localized — a typical buckling event involves only two to 
five particles. This is reminiscent of the Tl process found 
in a two-dimensional foam [33j. 
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FIG. 9: (Color online) Structural signature in higher order 
peaks of the pair correlation function, (a) Height of the second 
peak of pair correlation function, ^2, as a function of (j). (b) 
Height of the third peak of pair correlation function, ^3, as 
a function of 0. (c) Comparison of the amplitudes of 5^1 (0), 
^2(0) and ^3(0). The vertical dashed line indicates (/)/(/)c = 1- 



B. Displacement field and cluster structure 

More information can be obtained from the displace- 
ment field of the system. To visually illustrate the 
particle-displacement field, we subtract two images of the 
system at different times. Any stationary part of the sys- 
tem will appear black in such difference images since the 
individual pixels are identical in that region. However, 
if a particle moves, the image subtraction will produce 



8 



an area of crescent shape with positive values along the 
front boundary in the moving direction and leave a sim- 
ilar crescent area with negative values in the rear. Fur- 
thermore, we assign the negative values to zero (black). 
Hence, only front boundaries in the moving direction are 
indicated in the image. The curvature and area of the 
crescent show the direction and magnitude of a particle's 
displacement, respectively (Fig. fTOk-fTOtl). 

Fig. fTOk-HOli show the displacement field of a typical 
experiment. The displacement field is initially random 
(Fig. [TOk). As the system evolves toward higher pack- 
ing fractions, a coherent structure emerges (Fig. [TOb and 
[TOb). Particles tend to move outward around a few nuclei. 
In other words, a few clusters form in the system. Even- 
tually, when the system approaches the jamming point, 
one single cluster forms (Fig. flOli). 

How does the cluster structure emerge out of a random 
initial configuration? As mentioned above, since only a 
few pair of particles are in contact initially, the average 
displacement is small and the directions of displacement 
are random at the beginning (Fig. [TOk). As the size of 
particles increases, they begin to form local contact net- 
works — cluster structure emerges (Fig. [TOb). Encircled 
by its neighbors, a particle at the center of a cluster feels 
zero average force due to the balance of the compression 
from different directions. Hence, it stays in stationary 
and shows as a dark nucleus in the displacement field 
(Fig. [TOb). Meanwhile, the surrounding particles move 
away from the center radially. Due to inevitable initial 
density variations, the denser part of system will form 
clusters first. Since the displacement of particles at the 
edge of a cluster is linearly proportional to the size of the 
cluster (Fig.lHk and [8b), a larger cluster will grow faster, 
and therefore aggregate more particles. Thus, the clus- 
ters formed at early time in the denser part of system will 
quickly dominate the system and the small initial den- 
sity variation is amplified. When two clusters meet, the 
internal particles and especially those along the cluster 
boundaries rearrange and the two clusters merge into a 
single bigger entity (Fig. [TOb and [TOb). This merging of 
clusters continues until a single system-spanning cluster 
forms when the system approaches the jamming point 
(Fig.ttOH). 

The exact shape of evolving clusters depends on the 
initial condition of experiments. For the experiment 
shown in Fig. [TOl the initial density is higher in the cen- 
tral area of the system. Hence, the cluster structures 
appear first in that area and then expand outwardly. Al- 
though details depend on a packing's history, the exis- 
tence and development of cluster structures are robust. 
Unless over 10,000 particles in the system start out with, 
and maintain, exactly the same inter-particle spacings, 
which clearly is unrealistic, the system will always evolve 
into cluster structures at later time. Any small initial 
density variation will be amplified as the system ap- 
proaches the jamming point. 



C. Quantitative analysis 

To quantify the cluster structure, we measure the two- 
point correlation function of the displacement field, 

Cnnir) = ' ^ , (1) 

l^^(rl)-^(rl) 

i=l 

which is usually used to identify coherent structures in 
a system. Here, D{fi) is the displacement of the par- 
ticle i located at fi within a small time interval At, 
and Tij = — Tjl- The summations are over all the 
N particles in the system and the normalization factor 
^0 = ^i^j=i^{'^ij —t). Experimentally, the data are 
binned with a bin size of two thirds of a particle diameter. 
As one can see in Fig. [TTk . initially at low packing frac- 
tion (^/^c = 0.77) no correlation exists, Cgg{r > 0) = 0. 
As the packing fraction increases, clusters begin to form. 
Correspondingly, both the correlation and the correlation 
length increase fFig. [TTb). For example, the point where 
^DD crosses zero shifts to larger r as (/) increases. The 
correlation reaches a maximum when the mean square 
displacement {D'^{(j))) is largest. As {D'^ {(/))) falls near 
the jamming point (Fig. [8t), the correlation magnitude 
decreases but the length scale of the correlations is fixed 
at the size of the system (Fig. [TTb). One can of course 
choose other correlation functions for identifying cluster 
structures. The correlation function chosen here includes 
the information of both the direction and magnitude of 
displacements. We also measure the correlation of only 
the magnitudes of particle displacement, which is defined 
as 

1 ^ 

— ^ AI)(rl)AI)(r,)(5(r,,-r) 

CADAD{r)= ^ , (2) 

l^AI)(rl)AI)(rl) 

i=l 

where AD{fi) = D{fi) — (D) is the fiuctuation of the 
magnitude of particle displacement at r^. D{fi) is the 
magnitude of particle displacement, and {D) is the av- 
erage magnitude of particle displacement. We show 
Cadad{t) at different packing fractions in Fig. [121 As 
one can see, Cadad{t) shows qualitatively the same be- 
havior as Cgg{r). 

Cg^{r) (or Cadad{t)) shows a clear trend of the in- 
creasing of correlation between the displacements of par- 
ticles. It provides a good quantitative illustration of 
the emergence and the evolution of cluster structures. 
However, it is hard to extract the length scale of clus- 
ters from this function directly — Cgg{r) does not de- 
crease exponentially with r and there is no obvious fea- 
ture in C^^(r). Furthermore, due to the random noise 



FIG. 10: Cluster structure. Displacement field of a subsystem in the large cell (54.6 cm x 54.6 cm). The time interval for the 
displacement is At = 800 s. (a) = 0.65 at t = h, (b) = 0.74 at t = 1.30 h, (c) = 0.79 at t = 2.29 h, and (d) = 0.84 at 
t = 3.92 h. The white scale bar is 5 cm. 



in the displacement of individual particles, the correla- 
tion function shows a sharp jump from C^^(O) = 1 to 
Cgg{r = do). A particle always perfectly correlates with 
itself, but the correlation between different particles is re- 
duced due to random noise. Hence, there exists a sharp 
jump from the self-correlation at r = to the correlation 
of neighboring particles at r = (io- The jump is more se- 
vere when the signal/noise ratio is smaller at the begin- 
ning of the experiment or in the jammed phase (Fig. [TT]). 
(At the beginning, when particles do not touch and at 
the end when particles are jammed, the signal/noise ra- 
tio is ^ 0.) The random noise can be due to the non- 
uniform swelling of particles. If the shape of a parti- 
cle changes during swelling, the center of the particle 



may move slightly even without contact with other par- 
ticles. The particle- tracking algorithm can also induce 
some noise. However, that only happens when the system 
is deep inside the jammed phase where the inter-particle 
spacing is so small that it becomes hard to distinguish 
the boundary between two neighboring particles. 



Another way to quantify the cluster formation is to 
measure the projection of the relative displacement of 
two particles on the direction of their relative position 
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FIG. 11: (Color online) Correlation of the displacement 
of particles, C^^(r), before the mean square displacement 
reaches the maximum (a) and after the maximum near the 
jamming point (b). The arrows indicate the direction of in- 
creasing packing fraction. Values of (/)/(/)c are shown in the 
plots. The time interval for the displacement is At = 200 s. 
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FIG. 12: (Color online) Correlation of the magnitude of par- 
ticle displacement, CADAD{r), (a) before the mean square 
displacement reaches the maximum and (b) after the maxi- 
mum near the jamming point (b). The arrows indicate the 
direction of increasing packing fraction. Values of (/)/(/)c are 
shown in the plots. The time interval for the displacement is 
At 200 s. 



(Fig.im): 



1 ^ / " \ 

^'-W= AT S [{Din)-Dir-))-'^Uin,-r). 

(3) 

The bin size is again chosen as two thirds of a particle di- 



ameter. Dr{r) indicates on average whether a pair of par- 
ticles at separation r moves closer {Dr{r) < 0) or moves 
apart {Dj.{r) > 0). If there are no clusters and all the par- 
ticles move randomly, Dj.{r) will be zero. However, inside 
a cluster, any two particles cannot move closer (Fig.[T3k). 
Therefore, in a cluster with size ro, Dr{r < ro) > 0. 
When all particles are stationary, Dr{r) = 0. Figure [T3b 
shows Dr(r) at different packing fractions. As one can 
see, initially there is no coherent motion and Dr{r) is 
flat at zero. However, as the system evolves, Dr{r) be- 
gins to deviate from zero and a region with a positive 
Drir) clearly shows up implying the presence of clusters. 
Although the amplitude of this positive region shows a 
non-monotonic behavior, the position of the peak at / in- 
creases with (j). Thus, / can be used as a characteristic 
length scale of clusters. As shown in Fig. [T3h , / increases 
slowly from a few particles diameter at low (j) and in- 
creases rapidly to the size of the system when the system 
approaches the jamming point. The exact shape of /(0) 
depends on the initial configuration of the particles, but 
qualitatively all packs show the same behavior. 



D. Discussion 

In contrast to the length scales determined from dy- 
namic heterogeneities in supercooled liquids 0, IsHI . col- 
loidal suspensions |y, |36| and granular media |25M27|. 
the length measured here reflects the static structure of 
system. One might imagine that the average packing 
fraction inside a cluster is higher than that outside the 
cluster. The cluster structure observed here is due to 
the athermal nature of the system. Clearly, the clus- 
ter structure depends on the initial packing configura- 
tion. Although we prepare packs in a random way with 
unavoidable small density variations, the system always 
amplifies this initial variation into a cluster structure 
at later time. Such a history dependence and memory 



are typical of non-equilibrium systems |37l-l39j. On the 
contrary, in any equilibrium system without attraction, 
the initial density fluctuations will quickly be smeared 
out; unless the system approaches a second order phase 
transition, the correlations of the density fluctuations in 
equilibrium systems will remain small. The existence 
of a divergent static length scale in the glass or jam- 
ming transitions would be a hallmark for an underlying 
phase transition. However, until now, no unambiguous 
static length scale has yet been observed in any equilib- 
rium system near the glass or jamming transition. The 
present experiment suggests that at the jamming tran- 
sition a non- equilibrium system can produce a divergent 
static length scale. It should be noted that previous sim- 
ulation works on the jamming transition at zero temper- 
ature acquire the unjammed configuration by quenching 
the system from T = oo to T = Therefore, 
even though the system is at T = 0, the static config- 
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uration of the system is essentially the same as that at 
T = oo. Contacts between particles do not exist before 
jamming and static length scale cannot be found in these 
simulations. 



V. EFFECT OF FRICTION: MULTIPLE 
JAMMING POINTS 

As shown in the above section, the initial density vari- 
ation, no matter how small it is, will be amplified by the 
system en route to the jammed phase and the system will 
spontaneously organize into cluster structure. In this sec- 
tion, we shall investigate how this density variation influ- 
ences the signature of the jamming transition. We shall 
show that for a system with sufficient large initial den- 
sity variation, i.e. for a highly inhomogeneous system, 
friction plays an important role and multiple jamming 
points exist. 

Practically, the initial density variation can be effec- 
tively controlled by the total number of particles in the 
system. If enough particles are added at the beginning, 
the initial packing fraction of system will be high and 
the density variation will be small. In this case, although 
the cluster structure still emerges while the system ap- 
proaches the jamming point, the system will eventually 
show a clear jamming transition and therefore an unam- 
biguous structural signature of the transition as shown 
in Sec. III. However, if the sample is prepared at a low 
initial packing fraction, due to random vibrations dur- 
ing the sample preparation, there is a good chance that 
one part of the system is denser than the rest. The ini- 
tial density variation of the system will be much higher. 
Hence, we can easily prepare a highly inhomogeneous sys- 
tem by simply reducing the initial packing fraction. With 
a highly inhomogeneous system at hand, we want to ask 
how such a system goes through the jamming transition. 
Is the jamming signature of a highly inhomogeneous sys- 
tem the same as that of a homogeneous system? 

Figure [T4k shows the height of the first peaks of the 
pair correlation function, ^i, as a function of (j) for a sys- 
tem with low initial packing fraction (typically ^initial < 
0.60). Instead of a single pronounced peak, there are two 
major peaks located at ^ 0.74 and ^2 — 0.84 (and 
possible several smaller peaks). We can compare gi{(j)) 
with the mechanical and kinematic criteria of the jam- 
ming transition separately. As shown in Fig. [T4k . the 
pressure on the boundary begins to deviate from zero 
near the first peak ^1. However, if one compares gi{(j)) 
with the mean square displacement, (1^^), one finds that 
{D^) drops to zero at the second peak (j)2 (Fig. fT4b). 
Therefore, the system jams at (pi according to the me- 
chanical criterion for the jamming transition and jams at 
^2 according to the kinematic criterion. Why do the two 
well-defined criteria for the jamming transition occur at 
two different packing fractions? 
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FIG. 13: (Color online) Projection of the relative displace- 
ment of two particles on the direction of their relative position, 
Dr, and length scale, /, extracted from the cluster structure, 
(a) Illustration of the definition of Dr. A pair of particles with 
distance r (shown with red arrows) always move apart inside 
a cluster, (b) Dr as a function of r for different packing frac- 
tions. The time interval for the displacement is At = 200 s. 
Values of (/)/(/)c are shown in the plot. The horizontal dashed 
line indicates Dr = 0. (c) Length scale, /, extracted from Dr, 
as a function of ( 
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FIG. 14: (Color online) Height of the first peak of the pair 
correlation function, ^i, as a function of (j) for a system with 
low initial packing fraction, and comparison of ^i(0) with (a) 
the force along the boundary and with (b) the mean 

square displacement of particles gi{(j)) is shown on 

the left with black squares. The vertical lines mark the po- 
sitions of two major peaks, 01 and 02. (i) and (ii) indicates 
the positions where we show the displacement field in Fig. 1151 
F{(j)) is shown on the right of (a) with red line. The horizontal 
dashed line indicates zero force. (1)^(0)) is shown on the right 
of (b) with red circles. The time interval for the displacement 
is At = 800 s. 



To understand this, we plot the displacement field at 
two different packing fractions in Fig. [151 (i) at the first 
peak 01 (Fig. [T4b), and (ii) at the valley between 0i 
and 02 (Fig. [Hk). As one can see, at 0i the system 
already forms clusters (Fig. [T5b). From the local mean 
displacement (red arrows), one can see that the two clus- 
ters span the system from left to right. The cluster on 
the right pushes onto the force sensor at the boundary 
of cell. Clearly, the system jams and a force chain forms 
along this direction from left to right. However, there is 
still empty space in the lower right corner of cell. Af- 
ter the two clusters merge, the particles move together 
towards the lower right corner as shown in the displace- 
ment field near the valley (Fig. [T5b). The system jams 
globally at 02, which results in the dramatic drop of the 



mean square displacement. Therefore, the reason why 
there are two peaks in ^i(0) is that the system jams in 
two steps: it first jams locally along certain direction 
(from left to right) and then it jams globally. For any 
system like this, depending on where the force sensor is 
located, the pressure measured may deviate from zero at 
different packing fractions. But it should always happen 
before the global jamming point (02), as confirmed by 
the experiments. 

With frictionless particles, the system can only jam as 
a whole due to force balance throughout the entire sam- 
ple. Therefore, the multiple jamming points found in 
the system have to be due to friction between particles 
and between the particles and the boundary of cell. A 
frictional granular system can form very inhomogeneous 
structures such as force chains along one specific direc- 
tion and jam in that direction, but still have empty space 
in other directions. As the size of particles increases fur- 
ther, force chains will buckle under the increased stress, 
displacing particles into the less dense regions adjacent 
to the chains, and causing the system to unjam. As a 
result, the system approaches the transition to jamming 
in a series of steps. This picture is reminiscent of the sce- 
nario proposed based on theoretical considerations that 
rigidity emerges by successive buckling of force chains in 
glasses and granular matter 4o|, 41 1 . It should be empha- 
sized that frictional contact is an essential ingredient for 
the existence of force chains before the global jamming 
point. However, deep inside the jammed phase, where 
the rigidity is already well established, force chains can 
sustain without friction. In conclusion, with friction the 
picture of a single jamming transition at T = (Fig. [1]) 
has to be modified. Multiple jamming points may ex- 
ist in a frictional system when a highly inhomogeneous 
structure is present. 

To further test the above picture, we performed the ex- 
periment with small vertical vibrations applied to the sys- 
tem (Fig. [3]). By vibrating the system, any force chains 
formed before jamming are destroyed by relative slip be- 
tween particles. The vibration also helps one to demobi- 
lize the frictional contact between particles. If our picture 
of the relation between friction and jamming transition 
is correct, then by adding vibration the system should 
jam in one single step. In the experiment, the vibration 
is generated by a mechanical shaker in a tapping mode. 
Each tap is excited by one full period of a sinusoidal 
wave with the frequency = 30 Hz. The peak-to-peak 
acceleration of vibration is Tp-p = 1.51^, measured by 
an accelerometer attached to the cell. Here, g is the 
gravitational acceleration. The amplitude of vibration 
A = Tp-p/ (2uj'^) =0.84 mm is much smaller than the di- 
ameter of particles. Hence, particles only vibrate locally 
around their mean position. We tapped the cell once ev- 
ery 20 s before taking an image. A time interval of 2 s 
is allowed between shaking and taking an image, so the 
system is stationary when the image is taken. Six differ- 
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FIG. 15: (Color online) Displacement field of an inhomogeneous system in the smaller cell (16.5 cm x 16.5 cm), (a) Displacement 
field near the first peak of ^i((/)) at point (i) shown in Fig. [M] (b) Displacement field at the valley between the first and the 
second peaks at point (ii) shown in Fig. 1141 The time interval for the displacement is At = 800 s. Red arrows show the local 
mean displacement. The blue box on the right indicates the force sensor at the boundary of the cell. The white scale bar is 3 
cm. 



ent experiments all with low initial packing fractions are 
performed. As expected, none of the experiments shows 
the multiple- step jamming. All the systems show a single 
peak in gi{(j)). We show a typical result in Fig.[T6l Differ- 
ent from the results without vibration (Fig.[8h), the mean 
square displacement of particles, {D'^ {(/))) ^ decays mono- 
tonically (Fig. [T6|) . At low packing fraction, particles 
have more free room to vibrate and therefore have larger 
amplitude of displacement. The displacement amplitude 
decreases as the packing fraction of system increases and 
goes to zero when the system jams. Furthermore, the sys- 
tem under vibration does not develop any cluster struc- 
ture as it approaches the jamming point. As shown in 
Fig. [13 the correlation of the displacement of particles, 
C^^(r), keeps roughly the same shape as ^ ^c, which 
is clearly different from the non- vibrating case. When 
the cluster structure develops, C^^{r) increases signif- 
icantly (Fig. [TT]) . This confirms the argument that the 
cluster structure, and therefore the static length scale of 
the jamming transition, is due to the athermal nature of 
the system. Vibration thermalizes the system in a certain 
sense, and therefore destroys the cluster structure. 

VI. CONCLUSIONS 
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FIG. 16: (Color online) Structural signature of the jamming 
transition of a system subjected to small amplitude vibra- 
tions. The height of the first peak of pair correlation function 
gi is shown on the left with black squares. The mean square 
displacement of particles (D^((/))) is shown on the right with 
red circles. The time interval for the displacement is At = 200 
s. (f) is normalized by the packing fraction at the jamming 
point (j)c — 0.83, which is indicated by the vertical dashed 
line. 



In this paper, we systematically investigate the jam- 
ming transition in a 2D granular system. We show that 
there is a clear structural signature of the jamming tran- 
sition at zero temperature. The heights of the first and 
second peaks of the pair correlation function, ^i((/)) and 
^2(^)7 both show a maximum as the system crosses the 



jamming point. By measuring the pressure along the 
boundary and the displacement of particles, we show 
that this maximum coincides with the mechanical and 
kinematic criteria of jamming. Therefore, the structural 
signature we found here is a signature of the jamming 
transition. Although the amplitude of the peak does not 
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FIG. 17: (Color online) Correlation of the displacement of 
particles, C^^{r), for a system under vibration. Values of 



are shown in the plots. Dashed horizontal line indicates 
zero correlation. 



diverge due to the polydispersity of our particles, our 
experiment corroborates the results predicted in simula- 
tions with ideal frictionless particles at zero temperature 
p. This signature has already been used as a new cri- 
terion of jamming at finite temperature, where the me- 
chanical and kinematic criteria of jamming are hard to 
measure directly fi§\. Here, our measurement that shows 
the coincidence of the maximum with the onset of the 
rigidity provides an experimental basis for the criterion. 
Of particular importance is that this structural signature 
exists for real granular systems with frictional contacts. 
The structural signature found in experiments reflects 
the underlying singularity of the jamming transition. 

Friction may cause systems to jam in a series of steps. 
We found that if the initial packing configuration is 
highly inhomogenous, the system may jam in a certain 
direction but not in others. This phenomenon is directly 
related to the well-known force chain structure of a gran- 
ular system [5] . Our observation provides more insight on 
the relation between the heterogeneous force-chain struc- 
ture and the jamming transition in the presence of fric- 
tion We speculate that the jamming transition of 
frictional system is obtained through a continuous buck- 
ling of force chains in the different directions. 

It is also useful to note that ^i(0) is a very sensitive 
probe for the jamming of system. Even partial jamming 
can induce a peak in ((/>). By contrast with the kine- 
matic criterion of jamming, the average displacement or 
velocity of particles drops to zero only at the final global 
jamming point. For the mechanical criterion of jamming, 
depending on where one measures the pressure along the 
boundary, it may show jamming at different packing frac- 



tions. Therefore, the peak of which can be called 

the geometrical criterion, is a better jamming criterion 
for a system with friction. 

Due to its athermal nature, this system has a depen- 
dence on its initial particle configuration. It amplifies 
any small initial density variation and self-organizes into 
clusters. A static length scale extracted from this clus- 
ter structure reaches the system size when the system 
approaches the jamming point. Hence, we show a diver- 
gent static length scale in this nonequilibrium jamming 
system. 
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